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Motivation and Objectives 


¢ Motivation: 


—The probability of collision (Pc) between two Earth- 
orbiting satellites requires estimates of their orbital 
trajectories and associated uncertainties 


—Pc estimation requires processing conjunction state 
vectors and covariance matrices 


—CARA sometimes encounters non-positive definite 
(NPD) covariances that can potentially prevent or 
adversely affect Pc estimation 

¢ Objectives: 
—Investigate the frequency of NPD covariances 


—Implement method(s) to remediate NPD covariances 
when necessary 
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Geometrical Visualization 
of a 3x3 Covariance Matrix 


The 1-sigma ellipsoid of 
the PDF defined by the 
3x3 covariance matrix 


Major Eigenvector = V,_., 
Major Eigenvalue = 1,,., 
Semi-major axis = 6,4, = (Amax) 


Minor Eigenvector = V,,,, 
Semi-minor axis = o,,,, 


Medium Eigenvector = V,.4 
Ssemi-principal axis = One 
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Transition from Positive Definite 
to Semi-Positive Definite 3x3 Covariances 


Xd min = Amea Amin =Amea! 4 = Amin = Amea/ 100 Amin = 9 
Omin = Omed Omin = Omea/ 2 Omin = Omea/ 10 Omin = 9 


Positive definite (PD) covariances have positive o,,,, values 
Positive semi-definite (PSD) covariances have Zero o,,,, values 
Non-positive definite (NPD) covariances have imaginary o,,,, values 
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The Frequency of NPD State 
Covariance Matrices Over Time 


¢ Analysis of 830,509 events --- 6x6 ECI state covariance matrices at TCA 
— 2015-04-01 to 2017-04-06 ---- 6x6 ECI state correlation matrices at TCA 


¢ NPD covariances for primary 
objects decreased markedly 
just before 2016-05-01 
— Coincides with an increase in 
the number of significant 
figures used for covariances 


— Better precision leads to fewer 
NPD state covariances 


NPD Fraction (% ) 
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¢ This same pattern is seen for 
3x3 position covariances -600 400 -200 0 


— But at frequencies reduced by Time (days from 2017-01-01) 
a factor of 100 or more 


Omicron Doyle Hall | 7 


High Eccentricity vs. Low Eccentricity 
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¢ Analysis of 428,589 conjunctions 
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—High precision covariances 
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Covariance Requirements for Numerical 
Stability in Pc-Related Estimations 


¢ Mahalanobis distance estimation 
— The marginalized 3x3 relative position covariance A(t) = A,(t)+ A,(t) 
needs to be positive definite 
¢2D Pc estimation 
— A marginalized 2x2 projection of A(t.) = A,(t.4)+ A,(t.,) needs to be 
positive definite 
¢3D Pc estimation 
— The marginalized 3x3 relative position covariances A(t,) need to be 
positive definite at all ephemeris times t, used in the calculation 
¢ Monte Carlo Pc estimation 
— The full NxN state covariances at the sampling epoch time C(t,,) 
need to be at least positive semi-definite 


Pc-related calculations don’t always require fully PD state 
covariances for both objects, because they often use combined 
and marginalized covariances with reduced dimensions 
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Add an offset to all 
Description eigenvalues, and then 
reconstruct covariances 
using original eigenvectors 


-Relatively simple to 
Advantages implement 


-Assumes original 
Dis- eigenvectors can be used 
for matrix reconstruction 
hae -Offset applied to all 
eigenvalues, even if they 
don’t need remediation 
-Not constrained using any 
physics-based 
considerations 


Find closest PSD 
covariance and/or 
correlation matrix in terms 
of Frobenius norm 


-Mathematically well 
defined, as used by the 
financial industry 
-Algorithms and codes 
posted on-line 


-Most complicated to 
implement in software 
-Only designed to produce 
PSD matrices 
-The “closest Frobenius 
norm” criterion, while well 
mathematically defined, is 
not a physics-based 
criterion 


Three Methods Studied for Remediating 
Non-Positive Definite Covariance Matrices 


Clip the eigenvalues at a 

minimum limit, and then 

reconstruct covariances 
using original eigenvectors 


-Simplest to implement 
-Constrained by physics- 
based considerations 
- Produces fully PD 
position covariances 


-Assumes original 
eigenvectors can be used 
for matrix reconstruction 


oO _ “N.J. Higham, “Computing the Nearest Correlation Matrix—A Problem 
MICRON tom Finance,” IMA Journal of Numerical Analysis 22, 329-343, 2002. Doyle Hall | 11 


Pc Estimates for Conjunctions 
with Remediated NPD Covariances 


25338 _conj_32035_ 20161116 093255 20161112 231203 
TCA 6x6 ECI covariance status: C,=PD C,=NPD C,+C,=PD 


No covariance remediation 1.1137 x 107 
Spectrum shifting remediation 1.1137 x 107 
Higham remediation 1.1137 x 107 
Eigenvalue clipping 1.1137 x 107 


Most often, remediating 6x6 ECI covariances changes Pc negligibly 
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Pc Estimates for Conjunctions 
with Remediated NPD Covariances 


29499 conj_38481_20160727 045045 20160815 125105* 
TCA 6x6 ECI covariance status: C,=NPD C,=NPD C,+C,=NPD 


No covariance remediation 1.3217 x 10° 
Spectrum shifting remediation 1.3217 x 10° 
Higham remediation 1.3464 x 10° 
Eigenvalue clipping 1.3217 x 10° 


This ~1.8% Pc difference was the largest seen among 430,000 events 
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Remediating NPD Position Covariances 
using the Eigenvalue Clipping Method 


tT T TI] faz 0 O] fe vi > 
Figen -decomposition: C=/V, V, V, O A, O <« Vv > 
, + vi {0 0 Alle VO 


=() Remediation to PSD 
Clippedeigenvalues: 4...) =max[/,,4,,. | A.,. 20= —_ 

: 4 . >0Q Remediation to PD 
tT TF fF Ayn 0 0 ~~ ¥% 3 


Remediated matrix: CC, =|V, V, V, QO Axiom 9 < Vv, -> 
+ + ¥ 0 0 Aim)! |< Vv; > 


What is a sensible value for 1,,,, - the eigenvalue clipping limit? 
Can physics be used to constrain this value? 
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Collision Probability as a Function of o,,,, 


The collision probability represents an integral of the relative position 
PDF over the volume carved out along the path of the collision sphere 


Onin = 2 HBR# o,,, = HBR/10 Omin << HBR 


These three conjunctions produce similar Pc values 
(because they’ve carved out similar fractions of the PDF) 
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Collision Probability Visualization: 
Pc as a Function of o,,,, 


The collision probability represents an integral of the relative position 
PDF over the volume carved out along the path of the collision sphere 


The two conjunctions on the right will produce similar Pc values 
The one on the left will produce a smaller Pc value 
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Collision Probability Visualization: 
Pc as a Function of o,,,, 


The collision probability represents an integral of the relative position 
PDF over the volume carved out along the path of the collision sphere 


Pc values are insensitive to the o,,,, value whenever 0 < o,,,, << HBR 
This can be used to set a sensible eigenvalue clipping level 
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Conclusions 


¢ The frequency of NPD state covariance matrices decreased 
markedly in mid-2016 


— After increasing the covariance precision 
¢ NPD state covariances occur much more frequently for objects in 
high-eccentricity orbits 
—Likely due to covariance interpolation inaccuracies 


¢ Estimating Pc values doesn’t always require fully positive definite 
state covariances 


—Because the calculations use marginalized covariances 

¢ Remediating ECI state covariances doesn’t change Pc values 
—Three remediation methods produce equivalent results 

¢ Eigenvalue clipping can be used for remediation 


—It is the simplest to implement, and can be constrained when 
required using physics-based considerations 
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